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Ph. Abstract 



The systematic errors inherent in the QCD sum rule approach to meson 
wavefunctions are examined in the context of QCD in 1+1 spacetime dimen- 
sions in the large N limit where the theory is exactly solvable. It is shown that 
the truncation of high momentum modes induced in a lattice discretization au- 
tomatically produces a Chernyak-Zhitnitsky ^ type meson wavefunction. Such 
a truncation alters the balance of leading and higher twist terms in correlators. 
We find that the reliable extraction of (a few) higher moments is possible pro- 
vided a reasonably accurate uniform approximation to the Euclidean correlator 
over a suitable range is available, but that the extracted values are par- 
ticularly sensitive to the balance of lower and higher twist contributions. Un- 
derestimates of lower twist contributions or overestimates of the highest twist 
term may lead to too high values for the second and fourth moments of the pion 
wavefunction, suggesting a doubly peaked structure of the Chernyak-Zhitnitsky 
type. 
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1 Introduction 



Although lattice gauge theory remains the principal technique for extracting reliable 
nonperturbative information in 4 dimensional quantum chromodynamics (QCD), 
the method of QCD sum rules , in which physical results are extracted by combin- 
ing perturbation theory with nonperturbative information encapsulated in a small 
number of hadronic matrix elements ("condensates"), has also enjoyed considerable 
popularity. One of the more startling (and certainly unanticipated) results obtained 
by the sum rules technique is the prediction of Chernyhak and Zhitnitsky that the 
pion wavefunction in light-cone gauge has a doubly peaked, non-convex structure 
in which one of the quarks is most likely to be carrying the preponderance of the 
meson momentum. Some early lattice calculations found even larger values for 
the second moment, restoring the maximum probability to equally shared momenta, 
but producing negative values for the ground state wavefunction. Although a double 
peak wavefunction gives a better fit to the experimental data on the pion form factor 
at accessible energies, such an Ansatz has come under increasing scrutiny Q, Q, as 
it is unclear whether higher twist contributions (i.e. higher Fock states than simply 
a valence quark-antiquark) are really negligible in the accessible range, in par- 
ticular in the end-point regions which are necessarily emphasized by a wavefunction 
of Chernyak-Zhitnitsky type. 

In this paper we report the results of a careful study of the systematic errors 
intrinsic to the sum rules approach in an essentially solvable model with close resem- 
blance to 4 dimensional QCD, namely, the large N (= number of colors) limit of two- 
dimensional quantum chromodynamics. Quantized gauge theories become particu- 
larly simple in two space-time dimensions, a feature first exploited by SchwingerQ 
in his seminal paper on massless quantum electrodynamics. The analytic tractabil- 
ity of the nonabelian theory with gauge group SU(N) in the limit of large N was 
utilized by 't HooftQ in his study of the meson spectrum in 2-dimensional QCD. 
It turns out that this model allows us to study in great detail the accuracy of the 
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sum rule method and its sensitivity to various types of systematic errors in the in- 
put data in a situation where the correlators being studied and the spectrum and 
physical matrix elements which are essential ingredients of the procedure are under 
complete analytic or numerical control. 

One clue to the type of systematic error which can wreak havoc in a sum rules 
approach is provided by the observation that the CZ wavefunction appears auto- 
matically in a rather simple truncation of the theory: namely, in the boost of the 
lattice discretized Coulomb gauge wavefunction to light-cone gauge. In Section 2 
we explain how this comes about by showing how to connect a relativistic quark 
hamiltonian appropriate for an equal-time formalism with the 't Hooft Hamilto- 
nian appropriate for light-cone gauge. In Section 3 we review and extend some 
well-known properties of 2-dimensional QCD at large N. Section 4 contains some 1 
and 2-loop perturbative calculations which allow us to compute the required higher 
twist terms in the asymptotic expansion of the correlators studied in the sum rules 
approach to wavefunctions. In Section 5 we show how the method can be applied 
successfully, provided a reasonably accurate uniform fit to the Euclidean correlator 
is available over a suitable range. Precise control of logarithmic terms in the 
higher twist parts is found not to be crucial (although certainly more accurate re- 
sults are obtained if we include these terms correctly). In Section 6, the method is 
shown to fail seriously if the balance of lower and higher twist terms is altered (as it 
is in the lattice discretized situation discussed in Section 2). This may happen sim- 
ply because the higher dimension condensates depend on a high power of the QCD 
scale, which is not known with precision, or if unexpected nonperturbative terms 
(such as the infamous first infrared renormalon should happen to be present. 
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2 Light-Cone Wavefunctions in Lattice discretized 
QCD 

The bound state equation in 2 dimensional QCD in Coulomb gauge and in the 
limit Nc ^ oo has been thoroughly investigated by Bars and Green [^. Up to a 
unitary transformation the full Bethe-Salpeter wavefunction can be separated into 
two contributions 

, , (1 + 7o) ^ , (1 - 7o) ... 

(t> = (f>+ ^ l5 + (t>- ^ 75 (1) 

corresponding to a qq pair moving forward and backward in time. and (j)^ are 
defined by a coupled set of integral equations, 



P / 7 r^[Cip,k,r)<l)+ir,k) + Sip,k,r)^^ir,k)] (2) 



[E{p) + E{r — p) — r^](f>+{r,p) 

g^N r dk 

2tt J (p-k)^ 
[E{p) + E{r -p) + r°](/)_(r,p) = 

1 ■^^:r^2iC{p,k,r^.{r,k) + S{p,k,r)Mr,k)] , (3) 



where 



C{p,k,r) = cos^[G(p)-e(A:)]cos^[e(r-p)-e(r-A:)] , (4) 

S{p,k,r) = sin^[e(p)-e(A:)]sin^[e(r-A:)-e(r-p)] . (5) 

The functions E(p) and 0(p) are parametrizations of the quark self energy 

T,{p) = [E{p) cos @{p) — m] + j^[E{p) sin Q{p) — p] (6) 

and are determined by the self-consistent quark self energy equation. One can easily 
solve for E{p) and @{p) in the nonrelativistic (small coupling) limit, where eqs.(||) 
and (^) decouple and yield the ordinary non-relativistic Schrodinger equation. Like- 
wise one can solve for E(p) and Q(p) in the large momentum limit, finding 

p 



vr m 1 , 

2 p \p2 . 



(7) 



E{p) = \p\-W^ + i^\+0{-) . 8 
\p\ Z\p\ 



Bars and Green pointed out that in the boosted hmit r ^ oo and xp = r the 
coupled equations decouple as well, yielding a vanishing 0_ and one remaining 
integral equation for (/>_(_ similar to 't Hooft's integral equation Q obtained in light- 
cone gauge. In [10| it was shown how a Green's function in Coulomb gauge can be 



analytically transformed to light-cone gauge. This involves boosting the momenta 
of the wavefunctions and performing appropriate transformations in Dirac space to 
obtain the correct 7 matrix structure. The transformation formula reads 

SA'^Coulomb{^~^r,A~^p)SX^ = '^(^){r,p) '^■^ ^ Lightcone{r,p) , (9) 

with the Lorentz transformation 

(coshfw) — sinh(w) . 
' ^ I (10) 
— sinh(w) cosh(w) 

and the corresponding spinor matrix S\. Note that 

/ 

A~'p 



_i cj^oo exp(tj) 



V2 



P- 
P- 



ill) 



and the transformation ^ produces the dependence of the light-cone Green's func- 
tion on light-cone variables as expected. All the factors of exp(a;) which diverge in 
the limit uj ^ oo cancel in the full transformation of the Bethe-Salpeter equation 
(||),(^),(^) to the 't Hooft light-cone equation according to (^. 

For our purposes we may introduce a single relativistic Schrodinger equation 
(relativistic quark model) which intuitively describes the binding of two quarks in a 
linear rising potential: 

H(j){x) = (^\Jp^ + + \J {r — pY + tn^ — ma — nii, + "^1^1^ 4^{x) = e-(t>{x) . (12) 

This equation interpolates properly between the nonrelativistic limit (for r = 0) 
and the boosted high momentum limit of the exact Bethe-Salpeter equation (^),(^). 
In either case 0_ vanishes while equation (^) for cj)j^ plays the central role (with 
C = l, 5 = 0). 



In (12) the prefactor g^N/ir from the potential is removed and all masses and 
momenta are measured in units of the coupling constant gy/lN/ir), likewise distances 
in units of {gy/{N/-K))-^ . 

In the following the relativistic quark model ( p!^ ) is treated by discretization 
on a lattice. Plotting the resulting wavefunction in terms of the variable x = p/r 
yields the approximate (and in the weak coupling limit, the exact) light-cone result. 
The boost is accounted for by using a large value of total momentum r, meaning 
r » ma,mb. 

We perform the discretization in coordinate space before going back to momen- 
tum space by Fast Fourier Transformation. Introducing a lattice spacing a, the 
physical box size L = Na and ip as the integer lattice momentum corresponding to 
p, the latticized equation reads 



. sin ^]\ml + J r- sin ^i^!-^) ' + (13) 



^. a N J " ^ \a N , 

)-, OD 
K + J^ E |«min(n,7V - n)e'^'^'^'-'^>l'' = e4>i, . (14) 
iq,n=0 

As a result of the discretization, the differential operator id/dx corresponds to the 
momentum space operator (2/a) sin^Trip/N) instead of to the naive lattice momen- 
tum 2iTip/L. This difference becomes important for physical momenta (in our case 
the total momentum r) approaching the UV cut-off 2/a. 

We choose lattice sizes of = 512, N = 1024 and N = 2048 to investigate the 
ground state wavefunction for equal quark masses m = 1 shown in Fig(l). The value 
of r as well as the physical box size L is kept fixed for each value of m, such that 
with decreasing number of points the physical lattice spacing a = L/N increases 
and the resulting UV momentum cut-off decreases. The result is rather startling: as 
the ultraviolet cutoff (inverse lattice spacing) is reduced the characteristic convex, 
singly-peaked form of the ground-state meson wavefunction found by 't Hooft is 
replaced by a double peaked structure with a minimum at x = 1/2 - exactly the 
shape suggested by the QCD sum rule method of Chernyak and Zhitnitsky ||l|. 



or indeed, by some early lattice results [^! The cutoff of quark high momentum 
modes alters completely the balance of leading and subleading terms in the large 
behavior of correlators of quark-antiquark currents, of course, and we shall see 
below that it is precisely such an alteration that can lead to large systematic errors 
in the extraction of moments of meson wavefunctions in the sum rule approach. 

3 Review of some results from 2-dimensional QCD 

In the limit where the number of colors N is taken large, with g^N /-k held fixed, 
QCD in 1 space- 1 time dimension exhibits a discrete spectrum of stable mesons of 
mass (^=1)2,--)) ^"^^ duality is exact in the sense that 2 point correlators are 
meromorphic functions of expressible as sums over resonance poles. As pointed 
out by Callan et al.jl^], the large (f' behavior arising from the asymptotic freedom 
of the theory is intimately related to the asymptotic behavior for large k (meson 
excitation level) of the meson masses //fc and decay constants fk- 

A convenient test case for examining the sensitivity of the sum rule approach to 
truncations of the short distance expansion is provided by the two point correlators 
of the tower of light-cone operators 

Sn{x) ^ -^^(xh,CD_ri;{x) (15) 

which reduce in light-cone gauge (j4_=0) to 

Sn{x) = ^^(x)75(V_)Xx) (16) 
V jV 

The standard sum rule approach to the pion wavefunction due to Chernyak and 
Zhitnitsky j^] begins with a short distance expansion for a two point correlator of 
currents (or densities) which couple to the meson in question. It turns out that the 
tower of operators based on a pseudoscalar density has the same leading behavior for 
large q"^ in 2 dimensions as the two-point correlator of axial currents in 4 dimensional 
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-ln(-f )+0(l/g2), q2^_^2^oo (18) 



QCD. Hence, we study 

Mn{q^) = ^ / d'xe^"- < 0|r{5„(x)So(0)|0 > -{q^ ^ -/x^) (17) 

2 -"2 
vr ■ /i^ 

This correlator may also be written (in the large N limit) as a sum over resonance 
poles 

A: odd ^ /^fc 

k odd ^ ^fc 

where the squared meson masses are eigenvalues of the 't Hooft Hamiltonian 

H<l^k{x) = —^Mx) + P [' M^l^l^dy = 4Mx) (20) 
x{l-x) Jo [x-yY 

Here the parameter 7 is the bare quark mass (which we take equal to the antiquark 
mass throughout) squared in units where is set to unity. The resonance residues 
are quadratic in the moments 

fnk = /' %f^</.fc(x)dx (21) 

Jo x{l — X) 

with the normalization 

nl 

k2 







mx)dx = l (22) 
As shown in [pj|], for any given moment 



Itt A 1 
/„fc-— (l + -^ + 0(-^ln(A;))), k^^ (23) 

which, together with the asymptotic behavior derived originally by 't Hooft 

/z^ ~7r2/c + 2(7- l)ln(A:) + 0(1/A:) (24) 

implies the required free field behavior of Mn{q^) at large Q^. 
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The meson wavefunctions (f)kix) are not known analytically, but may be obtained 
to high accuracy numerically by expanding in a finite basis: 



D 



4>k{x) = Ckox'^{l - x)^ + \/2 ^ Ckj sin(7r(2j - l)x) 



(25) 



The parameter /3 is related to the quark mass by 7r/3/tan(7r/3) = 1 — 7. The first 



basis function appearing on the right-hand-side of (25) is necessary to properly treat 
the nonanalytic behavior at the endpoints x=0,l. In the above basis the 't Hooft 
Hamiltonian becomes a (L'+l)x(L'-|-l) matrix which may be numerically diagonal- 
ized by standard techniques. Adequate accuracy (typically to 5 significant places at 
least) in all the quantities computed from these eigenfunctions was obtained by tak- 
ing D=120. The groundstate pseudoscalar meson ("pion") wavefuntions for a range 
of quarkmasses (m^ =0.135,0.534,1.0,2.60) are shown in Fig. 2. These wavefunctions 



are in all cases convex, with a single extremum at x 
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Figure 2: Ground state (pion) wavefunctions in 2D QCD 



Once the (pkix) have been computed, accurate values for the moments fnk are 
readily obtained by numerical integration. The nonanalytic singularities at a; = 
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0, 1 are best avoided by using the following exact identities, which follow from the 
integral equation ( |20| ) 

fok = — [ <Pk{x)dx 
7 "'0 

^ ;^(l + ^ + 0(pln(fc))), fc-oo (26) 

where the coefficient Aq, which is related to the higher twist behavior of Mo(g^), 
will be determined analytically below. The second moment is determined in terms 
of the zeroth as follows (we use the property (j)k{x) = (j)k{^ — for the pseudoscalar 
mesons) 

1 

/2fe = / —r-s Al-2xf(t)k{x)dx 

Jo x{l-x) 

1 

= 2 / (- - A + Ax)(t)k{x)dx 
Jo X 

= (1 'i)fok 

while may be related to /ofc and a nonsingular integral which may be readily 
computed to high accuracy once the wavefunctions are known: 

hk = fok - ^^\flj^ x^<pk{x)dx (28) 

In the sum rules approach, the first three moments {fok, f2k, and f^k) are estimated 
and then used to draw conclusions about the shape of the pion wavefunction. We 
shall therefore restrict our attention to these quantities below. 

The sub dominant terms in the asymptotic expansion for large k of fnk are re- 
lated to logarithmic higher twist terms in Mn{q'^). For example, it follows from the 
resonance sum representation of Mo((?^) that, as = —(^ +oo 

AM^) + ^-f^" In % + 0(1/Q') (29) 

The pure l/Q"^ however (without a logarithm) is not determined by the large k 
asymptotics. Indeed, dropping a finite number of initial terms in ( p!9| ) clearly alters 
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the coefficient of l/Q"^, while leaving the large k behavior unchanged. In the fol- 
lowing section we will compute the next to leading twist contributions to M„(g^) 
directly from 2 loop perturbation theory. This is possible as a consequence of the 
superrenormalizability of the theory. 

The sum rules technique involves, as mentioned above, an approximate determi- 
nation of the first three moments of the groundstate pseudoscalar . In 2-dimensional 
QCD, the moments fnk vanish by parity in the pseudoscalar sector generated by the 
basis (p5| ) unless n is even. Thus we will be computing the quantities /oi,/2i and 
/41. In the large N limit they are obtained essentially exactly by the numerical 
procedure outlined above. 

4 Large Q'^ behavior of Correlators 

In the large N limit of 2 dimensional QCD, a superrenormalizable theory, the leading 
behavior of the correlators M„(g-^) in the deep Euclidean regime —q^ = Q'^^oo 
is given by the single one-loop graph (Fig(3a)) in which the quark-antiquark pair 
propagates freely . Keeping quark mass dependent terms, the Feynman integral is 
readily performed for general n and one finds at large 

Mr-P(Q2) ^ln^ + i(n-l)^ln^ 

4, ^ 1 m2 q2 



+ {;i + -,„-i)a-2|:— ))^+o(^i„;i,) 

- iQ-^i^') (30) 



The two-loop graphs depicted in Fig(3b,3c) have a leading asymptotic behavior ~ 
corresponding to dimension 2 operators in an operator product expansion. The 
self-energy (Fig. 3b) and exchange (Fig. 3c) graphs are not well-defined individually 
in light-cone gauge until a regularization procedure has been given for the gluon 
propagator. We shall take the momentum space gluon propagator to be -P(p-) = 
Re( ^ ^.^ )^, where the infinitesimal e can only be sent to zero after combining graphs 
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(3b) and (3c). 




<3. 



+ 




Figure 3: One and two-loop graphs contributing to the correlators Mn{Q'^ 



A lengthy but straightforward evaluation of the two-loop Feynman integrals for 
the self-energy graph Fig(3b) yields (for n = 0) the asymptotic behavior 



M^\Q-') -(^L^)(^ln(6) + ^ln( 



2,g'^N.. 1 



1 



IT TT 



Q 



^^ + 0(i,ln(-^^^^ 



m 



(31) 



while the exchange graph yields 



+ 0(— ln(^)2) 



(32) 



As expected, the regularization dependence cancels between the two graphs. The 
discontinuity of this corrrelator is infrared-safe so the mass-singularities in de- 
pendent terms must also cancel, as they do, thereby removing the logarithmic terms 

completely. The power singularity in the exchange graph is not physical and is 
in fact removed by the overall subtraction at = /U^ needed to define the overall 
amplitude. The total 2-loop contribution to this moment is thus (now, and hence- 
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forth, using units where = 1) 

Mo^-i°°P(Q^):.-^^ + 0(^) (33) 

The higher moments (we shall only need n = 2,4) may be similarly computed and 
one finds in all cases a pure power dependence 

Mri°°P(Q2) ^l^n-2)^ + 0(ij) (34) 

Since we are now in possession of the full ^ In^Q"^) contribution (arising solely from 
the one loop graph Fig(3a)) the subdominant asymptotic coefficients An (such that 
fnk Of 1 + 4f large k) are now determined 

-(7- 1-^(^0 + ^n)) = -{n-lh 

^ An = ^(2(l-n)7-l) (35) 

Thus A2 = Aq — in agreement with the exact identity (p^). 

Combining one and two-loop contributions, the asymptotic behavior of MniQ"^), 
through next to leading terms, is thus (after the overall subtraction at = iJ?) 
given by 




where we remind the reader that the dimensionless variable 7 is simply the squared 
quark mass m? in the natural units where = 1. These analytic results serve as 
a useful check on the numerical extraction of higher twist terms which we use below 
to implement the sum rules technique for the groundstate pion of this theory. 
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5 Extraction of Moments by the Sum Rules Technique 

In the standard sum rules approach to meson wavefunctions Q , the absorptive part 
of a current-current correlator is modelled by a sum of resonances with known mass 
values but unknown residues (typically at most the lowest two resonances are used) 
at low values and by the perturbative QCD expression at large q^. The Borel 
transform of this Ansatz (which is a smooth function, not a distribution!) is then 
fitted to the corresponding transform of an operator product expansion (OPE) for 
the same correlator including higher twist terms over a suitable range of the Borel 
variable. In 4 dimensional QCD two sets of higher twist operators are included, 
so that terms falling off like are present explicitly (a potential 1/Q^ 

contribution is assumed absent, which is a source of some controversy, to which 



we return below [g, 12]). The hadronic matrix elements appearing in the OPE are 
estimated phenomenologically. We shall follow the same procedure in 2-dimensional 
QCD. The only difference is the inevitable presence of l/Q^ terms arising from twist 
2 operators, so that the inclusion of next to leading and next to next to leading terms 
corresponds in this case to 1/Q^ and behavior. 

Define s = q^ (> in the timelike regime) and the Borel transform of f{s) as 

with the contour C as indicated in Fig(4). Thus the Borel transforms of 
{In(-s), iln(-s), i, ^,etc} are { 1 , ^""^^^1"^% , - ^^p: , etc} . Here 7e is the Euler 
constant. Writing (in the spacelike region —q^ = > 0) 

Mjm = ^,n(^)-{i|,n(^) + || 



(Q'-^')} + 0(l/Q'^) (38) 



The corresponding Borel Transform is given by 
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Figure 4: Contour defining Borel Transform 



+ 



^" (1 



M^'' 6 
+ 0{1/M' 



+ (ln(-^) + 1 - 7,)2) + ^(ln(-^) + 1 - 7e) + 



En 

(39) 



On the other hand, sphtting M„(s) into a contribution from low-lying resonances 
(with squared masses below some cutoff Sn) and a high energy perturbative piece 



d{s — Sn) ■ f ln{—§), one may also write 



Mn{M 



2 

— ( 

TT 



2\ ^ lg-5„/M2 



fc=Af 



Pnfc = fnkfok 



(40) 



(41) 



The duality cut variable Sn is allowed to float and is determined in the fitting 
procedure when the right hand sides of (|39|) and (|40| ) are matched, as are the residues 
Pnk- We shall take N=2 and fix the masses of the two lowest pseudoscalars at their 
exact values obtained by diagonalizing the 't Hooft Hamiltonian (pO|). Information 
about the shape of the groundstate meson ("pion") wavefunction is obtained from 
the ratios ^^i = ^ and = ^. 

POl JOl POl JOl 
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We first examine the accuracy of tlie sum rules procedure on tlie assumption that 
the high momentum behavior of M„((5^) is known very accurately over a range of 

where a sizable fraction of the variation is due to the higher twist contributions. 
Specifically, the procedure involves the following steps: 



1. The resonance sum representation is used to calculate MniQ'^) to high accu- 
racy over a range fj? < < 200, where the subtraction point has been chosen 
throughout as ^u^ = 5 (in units where the squared mass scale intrinsic to the 
theory, namely is set to unity). For example, for the quark mass value 

= 0.534, with fi\ = 4.59, we take /i^ = 5. The resulting values are then 
fit over a range Qg < < 200 (there is almost no dependence on the upper 
cutoff, once it is chosen reasonably large) to the asymptotic expansion (^), 
allowing us to extract the coefficients ...Ey^. The low point of the fit 

range Qg is adjusted to obtain the best fit , and typically one finds —20- 
30, with a rms deviation < 10^^. As a check that this procedure is yielding 
sensible results, we can compare the coefficients An,Bn obtained from the fit 
with the analytic values obtained in the preceding section. The results for 
nig =0.534 are shown in the top of Table 1, for moments n =0,2,4, with the 
analytic (i.e. large N) values in parentheses. Evidently, this fitting procedure 
reproduces the twist 2 coefficients to about 10%. We shall argue below that 
the precise values of the coefficients are in fact not too important, as long as 
one has a good uniform fit to the correlators Mn over a large enough range. 

2. Once the coefficients A^, Bn, .., En are known the Borel transform (|39| ) can be 
computed over any given range M|f < < 200. Then a fit is performed 
by matching the Ansatz ( |39| ) to ( ]40| ) over a range of values where the 
higher twist contribution is appreciable. A convenient choice yielding perfectly 
reasonable results is to take Mg ~ ^■^(= 5). The fitting parameters in the 
Ansatz (^) are Sn, the duality cut, and the residues of the two lowest mesons 
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Quark mass (squared) = 0.534 



n (Moment) 




Bn 


Cn 




En 


fnl 





0.596(0.680) 


2.02(2.114) 


-0.398 


2.51 


1.98 


3.47(3.41) 


2 


-0.775(-0.680) 


-1.948(-2.201) 


3.303 


4.742 


1.269 


1.60(1.82) 


4 


-2.11(-2.040) 


-3.31(-3.796) 


3.81 


6.92 


4.56 


1.37(1.40) 


Quark mass (squared) = 2.60 





3.27(3.31) 


0.61(0.13) 


-1.84 


-13.7 


-2.89 


4.87(4.83) 


2 


-3.69(-3.31) 


-4.17(-5.48) 


-3.06 


1.93 


4.11 


1.75(1.58) 


4 


-9.58(-9.93) 


2.95(2.16) 


-47.4 


-18.9 


0.83 


1.61(0.97) 



Table 1: Asymptotic coefficients and moments (exact values in parentheses) 



Pni,Pn2,n = 0,2,4. From the residues pnk one easily solves for the moments 
fnk (cf (^T|)). The results for the moments are shown in the last column 
of Table 1, together with the exact values in parentheses. The first three 
moments are reproduced fairly well, to 10% accuracy. The same procedure 
applied to a heavier mass quark ("onium" type meson, with =2.6) leads to 
OPE coefficients and moments shown at the bottom of Table 1. The first two 
moments are still given reasonably well, but the n = 4 moment is too large. 

At this point it is appropriate to point out that accurate control of logarithmic 
terms in the asymptotics of M„((5^) is not essential for this procedure to give rea- 
sonably accurate results. Indeed, the Borel Transform (|38|) to (39) is evidently a 



linear mapping which becomes a bounded one (relative to the norm, which is 
relevant here as we perform least square fits throughout) once finite intervals in 
and M^, and a finite truncation of the OPE, are chosen. Thus the results obtained 
for Pnk are really quite insensitive to the precise values of the coefficients An, - -En, 
provided only that the asymptotic form is a good representation in the mean of the 
exact Mn{Q'^) over a suitable range. For example, if we use a pure power fit for 
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the higher twist terms 

M„(Q^).^,„,^)-(§ + § + | + f} (42) 

which gives a fit with rms deviation < 10~^ in the range < Q'^ < 200 (considerably 
worse than the fit obtained through twist 4, but including logarithmic terms), one 
finds /oi =3.38, /21 =1.95, and /41 =1.70. Only the last moment , n=4, seems to 
suffer (to the extent of about a 20% error) from the use of a pure power Ansatz 
for higher twist terms. Typically the QCD applications |jl| employ only pure power 
dependencies in higher twist terms. 

6 Truncation Sensitivity of the Method 

We saw previously (cf Section 2) that as simple a modification of the structure of the 
theory as the introduction of a lattice high-momentum cutoff can lead to a profound 
modification of the form of the ground-state meson wavefunctions. In the language 
of the Borel transform, the effect of such a cutoff is readily determined. A short 
calculation shows that the one-loop contribution to the Borel transform M„(M^) in 
the presence of a UV cutoff A = vr/a (a the lattice spacing) is given by 

f>7r/a 
t/q 

where 



M^'°°p(m2) = ^ r^'^ e-'^^'l'^'El-^iEl - m^yl^dk (43) 



i?fc^y^sin2(^)+m2 (44) 

In the continuum limit a ^ the integral ( ^3| ) becomes divergent for 00 and 

the large behavior is amplified to a constant term. For example, Mq{M'^) is 
given explicitly in the continuum limit by f ^e-^™'/*^' (i^o(i^) + Ki(^)) ^ 
;|, IVP 00. By contrast, in the presence of the UV lattice cutoff, the leading 
"twist-0" term is removed entirely, and the asymptotic behavior begins at order 
l/M^. So we can certainly expect trouble whenever the balance of lower and higher 
twist terms is altered. 
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We saw in the preceding section that the sum rules method appears to give 
reasonable results once a uniformly accurate fit to the M„ correlators is given. We 
shall now show directly that it can fail quite substantially if the balance between the 
various higher twist nonperturbative contributions is altered. In the applications of 
the method to 4-dimensional QCD the overall size of these higher twist terms 
is determined by a phenomenologically estimated condensate (i.e. expectation of 
a hadronic composite operator), while the moment dependence is obtained from 
a perturbatively computed coefficient function. Assuming the absence of the first 
infrared renormalon (on which more below) the twist 4 and 6 contributions are 
estimated from sum rules for charmonium [|l3| and from PCAC Thus different 
systematic errors are possible in different higher twist terms. Moreover, as these 
condensates involve the scale of the theory Aqcd to high powers, a relatively small 
ambiguity in the scale of the theory can alter the balance of higher twist terms 
quite substantially. We can examine the effect of a similar systematic error in the 
2-dimensional model by applying a scale factor to the twist 4 term, leaving the 
perturbative and twist 2 terms unchanged. Imitating the procedure used for QCD 
as closely as possible we first extract the best fit to the exact correlators M„((5^) (for 

=0.534) using pure powers only (through l/Q^), and then plot the sensitivity 
of the ratios /21 / /oi , /41 / /oi to an overall rescaling of the twist 4 contribution (see 
Fig. 5). Evidently, the higher moments increase steadily with the overall scale of the 
twist 4 term. In the QCD case, the considerably higher values found for the second 
and fourth moments (as compared to the asymptotic wavefunction oc x(l — x)) were 
interpreted by Chernyak and Zhitnitsky as evidence for a non-convex, doubly 
peaked wavefunction. 

It was recently argued |^ |l^ that there is no known rigorous argument to exclude 
the presence of an intrinsically nonperturbative 1/Q^ contribution to the coefficient 
function of the identity operator, and hence to the correlators M„(Q^) (this issue is 
intimately related with the location of the first infrared renormalon in the theory). 
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Figure 5: Sensitivity of moments to uniform rescaling of twist-4 terms 



In fact, if we return to the full logarithmic fits summarized in Table 1, and assume 
ignorance of the twist 2 terms by setting the coefficients An, Bn=0, we find perfectly 
consistent fits to the duality Ansatz (^0|). It should be emphasized that the fitting 
procedure does not provide any internal evidence that important contributions are 
missing. Indeed, the quality of the fits (using rms deviation as a figure of merit) 
actually improves by almost an order of magnitude! And the optimal fits find the 
duality cut variable (which is allowed to float in the fitting procedure) S„ settling 
at a self-consistent value, namely < Sn < A*! (since the lowest two mesons are 
included explicitly in the resonance sum part). The effect on the moments however 
is dramatic, and goes in the same direction as an increase in twist 4 contributions 
relative to lower twist terms discussed above. One now finds /21//01 — 1.0 and 
/41//01 — 1-1 implying a wavefunction with a minimum at x = 1/2, rather than the 
correct convex and singly peaked result. 

As regards the implications of our results for 4 dimensional QCD, it must be 
admitted that the prospects for a really firm nonperturbative resolution of the pion 
wavefunction seem rather dim at present. From the perspective of QCD sum rules, 
a reliable result would seem to depend on reasonably accurate uniform control over 



21 



the correlators over a wide range. First principles calculations starting from 
lattice QCD, on the other hand, also seem to introduce truncations of the theory 
which lead to potentially dangerous distortions of wave-function structure. Further 
calculations to study the systematic lattice errors in both correlators and moments 
on 4 dimensional (quenched) lattices of the size presently available are in progress. 
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